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ABSTRACT 

In this paper we attempted to present a policy to control an immune cell-tumor cell-normal cell-drug model 
proposed by Pillis et all. The drug administered to the patient in the form of chemotherapy is assumed to be time dependent 
and follows a definite rule. It is also assumed that the drug kills all types of cells. In this paper we assumed that the drug 
administration follows either of the three different mathematical laws viz. (1) Logistic law, (2) Exponential law and (3) 
Oscillatory law. Stability analyses of the tumor free critical points are done to find a range for the amount of drug to be 
administered to the patient. 
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INTRODUCTION 

The growth of cancerous tumor is a very complex phenomenon that includes many biological interactions. Till 
date many authors has proposed many tumor cell growth models and has suggested various control policies including 
treatment such as surgery, radiotherapy, drug therapy (Chemotherapy) [3], [4], [5], [6], [7], [8], [9], [10], [11], [13], [14]. In 2000, 
Pillis et all proposed a tumor growth model that involves immune cells, tumor cells and normal cells [12]. In the model 
proposed by Pillis et all, the immune cells and the tumor cells compete in a predator prey fashion whereas the normal cells 
and the tumor cells compete for available resources. The drug administered is assumed to kill all types of cells with 
different kill rates. In [2], the authors studied an almost similar model with the difference that the normal cells were not 
considered to interfere in the tumor growth model, with the drug administration following any of the three laws namely the 
logistic law, the exponential law and the oscillatory law. In this paper we have considered the tumor model suggested by 
Pillis et all [12], in which we have assumed that the tumor cells has no effect on the growth of normal cells as the number 
of tumor cells are negligible as compared to that of normal cells and so the normal cells does not need to compete with the 
tumor cells for the available resources whereas the normal cells effects the growth of tumor cells. This is possible as 
because we have considered the model for small sized tumor. We have assumed the drug administration in terms of 
chemotherapy follows one of the three basic mathematical laws viz. the logistic law, the exponential law and the oscillatory 
law. 

Main Study 

The model proposed by Pillis ET all in [12] is 

- = s + — - cJT - d-,1 -0,(1- e~ v )I 

dt a +T 1 1 IV J 
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f t = rj{ 1 - bj ) - c 2 /r - c 3 7W - a 2 (l - e-^r 
^ = r 2 JV(l -A/)- c 4 7W - a 3 (l - 


dl/ 

dt 


f(V) 


Where, /V (t) and V (t) denotes the number of immune cells, tumor cells, normal cells and the amount 

of drug administered at time t respectively. The units of cells are normalized by taking the carrying capacity of normal 
cells to be one. 


Here, s = the constant number of immune cells present in the body, 
er = Steepness coefficients. 

p = Recruitment rate of immune cells stimulated by the presence of tumor cells 
d 1= Natural death rate of immune cells. 


r 4 = Intrinsic tumor growth rate. 

= Tumor population carrying capacity. 

%, a 2 and a 3 are the kill rates of immune cells, tumor cells and normal cells respectively due to drug 
administration. 

In the above model we have considered a small sized tumor so that the number of tumor cells is negligible as 
compared to that of normal cells. So the tumor cells has no significant effect on the growth of normal cells. Thus we can 
choose the value of the parameter c 4 to be zero. 

As studied by the authors in [2], we also assumed three different drug administration model as given below: 

The drug equation as per logistic growth is as follows: 


%=a 3 V( 1-fcV) 


The drug equation as per exponential growth is as follows: 


— = a 3 e — y 3 
at 6 ' s 


The drug equation as per oscillatory growth is as follows: 


( 1 ) 


( 2 ) 


-=a 3 Sin(p 3 V) 


( 3 ) 


Where a 3 — intrinsic rate of drug application, 1 / fi 3 — maximum drug carrying capacity, a 3 = a constant rate of 
reduction of drug. 

Since the main objective of drug administration is to make the tumor size zero so we have considered the tumor 
free equilibrium points and analyzed their stability. 
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■ = 0 — > V = 0 or V 


- 1 


/ft 


dt ' 

When.F = 0, — = 0 gives T = 0 or T = ri ~ C2 '~ C3N 

dt b b ± r 4 


When, V = 0,T = 0,^=0 gives N = 0 or N = 1 


But N — 0 is the dead case and so it is discounted. 


dl 


So when, 1/ = 0, T — 0, N = 1, — — 0 gives I = s J ^ — I 0 (say ) 


Thus a tumor free equilibrium point is A(I 0 , 0,1,0), which is also drug free. 


Again for V — p , ^ = 0 gives T — 0 or T = 


1 //? 3 


r 1 -c 2 I-c 3 N-a 2 (l-e ^ 3 ) 

biri 


a 3 (l-e ^3) 


For V = l /p , T = 0,^=0 gives N = 0 or N = 1 - — — = N^say) 


When, V = 1 /p , T - 0, N = N ± then ^ = 0 gives / 


■ = h(say) 


d 1 +a 1 (l-e ^ 3 ) 

Thus another tumor free equilibrium point is B (I 1 , 0, /V-, , ^ J p ). 

The Jacobian matrix of the tumor growth model for the logistic drug administration at the equilibrium point A is, 


Ja 


( —d 1 ^ — c 4 I 0 0 O-iIq ^ 

0 A - c 2 /„ - c 3 0 0 

0 0 — r 2 -a 3 

\0 0 0 a 3 J 


The eigen values of J J; are: 


/li — d 1 , /l 2 — c 2 / 0 c 3 , /l 3 — r 2 and A 4 — ct 3 . 

The equilibrium point A is stable only when A 1 ,A 2 ,A 3 and A 4 are all less than zero. But A 4 = a 3 < 0 is not 
possible because negative drug administration is not a realistic case Thus there does not exist any tumor free drug free 
stable equilibrium point i.e. the tumor size can not be diminished to zero without drug administration and this is a very 
realistic case. 


The Jacobian matrix of the tumor growth model for the logistic drug administration at the equilibrium point B is, 


/ -d 1 -a 1 (1-e ^ 3 ) 


0 — the ^3 L 


J, 


L _ 
B 


i 1 


0 r i - c 2 /j - c 3 N, - a 2 (l - e ft 3 ) 0 

i0 o (1 -2^)^ -a 3 (l -e -a 3 e 

\° 0 0 -« 3 


7 


The eigen values of Jj] are: 

A x = -d 4 - %(! - e ^3), 
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A 2 = ri - C 2 h - C ? N 1 - a 2 (l - e ^3), 


h = (1 ~2N 1 )r 2 -a 3 (l - e ^ 3 ) 


and A 4 = —a 3 . 


Since A 4 = — a 3 < 0 is obvious so for B to be stable, we must have, A r < 0, A 2 < 0 and A 3 < 0. 

1 CL CL 1 XX 

A - 1 < 0 gives — > Loq( — — ) , which holds for any B 3 because — — < 1 => sq — < 0 but — > 0 ( — 

1 6 /? 3 aK a 1 + d 1 J J ai+di ft P3 


represents the maximum amount of drug administration and so it is positive). 


As < 0 gives /? 3 > 


- — — . But this happens only whena 3 — r 2 > 0. 


Again A 2 < 0 gives /? 3 < X_ — _ where m — r x — c 3 — — and n — — + a 2 

1 r 2 

Thus for the tumor free equilibrium point B to be stable, we must restrict the maximum amount of drug 
1 

administration — by the rule 

/?3 


[ < P 3 < , if a 3 - r 2 > 0 


Log{ —^3 — ) 


Log{ ) 

a K n-m J 


And 


0 <P 3 < 


Log ( ) 

a K n-m J 


if a 3 -r 2 < 0 


The following are the plots of number of immune cells, tumor cells and normal cells Vs. time respectively for the 
parameter values within the restricted range for the stability of the equilibrium point B. The plots clearly shows the 
stability of the equilibrium point B. The parameter values we chose to plot the graphs are s — 0.05, p = 1, a = 0.3, 
Cj = 0.2, d 4 — 0.2, a 4 — 0.2, r 4 = 0.4, b 4 — 1.5, c 2 = 0.3, c 3 = 0.2, a 2 — 0.5, r 2 = 0.35, c 4 — 0, a 3 — 0.25, a 3 = 0.5, 
/? 3 = 1.6 and the initial points are taken as 7(0) = , T(0) = 10~ 5 , 7V(0) = 0.9, V(0) = 10 -5 . 
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Figure 1 


Exponential Drug Administration Model 

37 = 0 => V = 7ft MS = V( s ay) 


c 1/ T/' dT n ■ r v c\ r p r 1 -c 2 I-c 3 N-a 2 (l-e v ) 

For V = V , — =0 gives T = 0 or T = 

dt b 


For V = V',T = 0,^=0 gives N = 0otN=1- 


N — 0 is the dead case and so it is discounted. 


When, V = V, T = 0, N = TV' then ^ = 0 gives / 


a 3 (i-e~ y ) 
r 2 


= N'(say ) 


-777 = /'(say) 


Thus the only tumor free equilibrium point is P(l ' , 0, N', I/'). 

The Jacobian matrix of the tumor growth model for the logistic drug administration at the equilibrium point P is. 


J L p = 


\ 


V , 

(1 — 2N')r 2 — a 3 (l — e~ v ') -a 3 e~ v N' 

-oc 3 P 3 e-W ) 


( -di-fliC l-e~ v ') —~cj 

0 r 4 — c 2 l' — c 3 N' — a 2 (l — e~ v ) 

0 0 

\0 0 0 

The eigen values of Jp are: 

A x — —d 1 — a 1 ( 1 — e ~ v '), 

A 2 =r t - c 2 l’ - c 3 N’ - a 2 ( 1 - e~ v '), 

A 3 = (1 - 2N 1 )r 2 - a 3 ( 1 - e“ v '') and A 4 = -a 3 /S 3 e~^ v '. 

Since A 4 = — a 3 /? 3 e _ ^ 3V ' < 0 always holds so for P to be stable we need only A 4 < 0, A 2 < 0 and A 3 < 0 


A x < 0 Gives /? 3 < 


f) 


L °9& f) 

A 3 < 0 Gives /? 3 > ; — , this holds only when a 3 — r 2 > 0 

L ° a ^a 3 -r 2 ^ 


(i) 


(ii) 
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Again X 2 < 0 Gives ft < -4— , where m — r x — c 3 and n = 1- a 2 (iii) 

1 r 2 

Thus for the tumor free equilibrium point P to be stable, we must restrict the maximum amount of drug 
1 

administration — by the rule 

P3 

Log( —) Log(—) Log(—) 

— 7 y < ft < Minimum{ — 7 V £ 3 / 7 , — 7 ^ 7 } if a 3 - r 2 > 0 

Msrsr) L °a{^) 

And 


0 < ft < Minimum} 


M ft) 


L °a{^) 

’ L °a{^) 


} if a 3 -r 2 < 0 


The following are the plots of number of immune cells, tumor cells and normal cells Vs. time respectively for the 
parameter values within the restricted range for the stability of the equilibrium point P. The plot clearly shows the stability 
of the equilibrium pointP. The parameter values we chose to plot the graphs are s — 0.05, p = 1, a — 0.3, c 1 = 0.2, 
ft = 0.2, a t — 0.2, r t — 0.4, b 1 = 1.5, c 2 — 0.3, c 3 = 0.2, a 2 = 0.5, r 2 = 0.35, c 4 = 0, a 3 = 0.25, a 3 — 0.5, y 3 = 0.05, 
ft = 3 and the initial points are taken as 7(0) = s / ^ , 7(0) = 10 -5 , A1(0) = 0.9, F(0) = 10 -5 . 




Figure 2 


Oscillatory Drug Administration Model 

kn / — 1 /* 


dV 

dt 


— 0 => V — Kn / p — V* (say}, where k is a positive integer 


t dT r 1 -c 2 I-c 3 N-a 2 (l-e v ) 

Fori/ — V , — =0 gives T — 0 or T = 

dt b b 


Fork = V*,T = 0,^ = 0 gives JV = 0orlV = l- a3{1 ^ } = N*(say ) 
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As N — 0 is the dead case and so it is not feasible and discounted. 

r] I c 

When, V = V*, T = 0, N = N* then — = 0 gives / = = I*(say) 

dt 6 d 1 +a 1 (l-e~ v ) v J J 

Thus the only tumor free equilibrium point for this model is Q(l* , 0, N* , V*). 

The Jacobian matrix of the tumor growth model for the logistic drug administration at the equilibrium point Q is. 


Jq 


L _ 


( d 4 6 

0 ry — c 2 I* - c 3 N 
0 


7 *) P -- Cl r 


a 2 (l 


0 —a 1 e~ v *r \ 

0 0 
(1 — 2N*)r 2 — a 3 (l — e~ v *) - a 3 e~ v *N * 

0 -a 3 p 3 Cos(P 3 V*) ) 


The Eigen values of Jq are: 

A 1 — —d 4 — a 4 {l — e~ v ), 

^2 = ri - c 2 r - c 3 N* - a 2 (l - e- 1 "), 

A 3 = (1 - 2JV*)r z - a 3 ( 1 - e~ v ') andA 4 = -a&CostfoV). 

For Q to be stable,^ < 0, A 2 < 0, A 3 < 0 and A 4 < 0 must be satisfied. 

A 4 < 0 => Cos(p 3 V *) < 0 => Cosjkn ) < 0 => k is odd 

, Log{ ~ j ■) ^ 

A, < 0 gives — > — — 1 , which holds for any /?, because — — < 1 => ^ — < 0 but — > 0. 

1 6 ft kn : a 1 +d 1 Log {- ft 

A 3 < 0 gives /? 3 > ; — kn a3 , this holds only when a 3 — r 2 > 0 

Loa ^a 3 -r z > 

* ■ i /-v • n fc?T , C?S t CnCL't 

Again a 2 < 0 gives p 3 < n — , where m = r x — c 3 and n = 1- a 2 

L°9( n -m) 1 72 


(b) 

(c) 


Thus for the tumor free equilibrium point Q to be stable, we must restrict the maximum amount of drug 

i 

administration — by the rule 

/?3 


kn 


kn 


L °g{- 


< ft < n ) if a 3 -r 2 > 0 


And 


0<ft < 


Log(——) 
a K n-m J 


if a 3 -r 2 < 0 


Here k is an odd positive integer. Since we are interested in minimum drug administration so that the normal cells 
gets less effected, so we chose A: = 1. Then V* — n / p and the relation between the parameters become 


L °g{- 


<ft < 




if a 3 -r 2 > 0 


And 


www.iaset.us 


editor @iaset. us 


50 


Ranu Paul 


°<ft <^pqi/«3-r 2 <0 

a K n-m J 

The following are the plots of number of immune cells, tumor cells and normal cells vs. time respectively for the 
parameter values within the restricted range for the stability of the equilibrium pointQ. The plot clearly shows the stability 
of the equilibrium pointQ. The parameter values we chose to plot the graphs are s = 0.05, p — 1, a — 0.3, c 1 — 0.2, 
d-i = 0.2, a-L = 0.2, r t = 0.4, b t — 1.5, c 2 = 0.3, c 3 = 0.2, a 2 = 0.5, r 2 = 0.35, c 4 = 0, a 3 = 0.25, a 3 = 0.5, /? 3 = 4.5 
and the initial points are taken as 7(0) = s , 7(0) = 10 -5 , A1(0) = 0.9, P(0) = 10~ 5 . 




CONCLUSIONS 

The stability of all the tumor free equilibrium points of the immune-tumor- normal cell model subjected to the 
three different drug administration laws are analyzed. It is found that for all the three models there does not exist any stable 
tumor free drug free equilibrium point whereas there exists tumor free stable equilibrium points satisfying certain 
conditions for the maximum drug dose. 
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